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Abstract 

The contact between a liquid and an elastic solid generates a stress vector depend- 
ing on the curvature tensor in each point of the separating surface. For nanometer 
values of the mean curvature and for suitable materials, the stress vector takes sig- 
nificant amplitude on the surface. Although the surface average action of the liquid 
on the solid is the hydrostatic pressure, the local strain generates torques tending 
to regularize the surface undulations and asperities. 
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1 Introduction 



In 1977, John Cahn gave simple illuminating arguments to describe the inter- 
action between solids and liquids. His model was based on a generalized van 
der Waals theory of fluids treated as attracting hard spheres [I]. It entailed as- 
signing an energy to the solid surface that is a functional of the liquid density 
at the surface. It was thoroughly examined in a review paper by de Gennes [2]. 
Three hypotheses are implicit in Cahn's picture: i) The liquid density is taken 
to be a smooth function of the distance from the solid surface, that surface is 
assumed to be flat on the scale of molecular sizes and the correlation length is 
assumed to be greater than intermoiecular distances; ii) The forces between 
solid and liquid are of short range with respect to intermoiecular distances; 
Hi) The liquid is considered in the framework of a mean- field theory. This 
means, in particular, that the free energy of the liquid is a classical so-called 



Email address: henri.gouin@univ-cezanne.fr (Henri Gouin). 



Preprint submitted to : Comptes Renclus Mecanique, vol. 337 (2009) 



gradient square functional. 

The point of view that the liquid in an interfacial region may be treated as 
bulk phase with a local free-energy density and an additional contribution 
arising from the nonuniformity which may be approximated by a gradient 
expansion truncated at the second order, is most likely to be successful and 
perhaps even quantitatively accurate near the liquid critical point [3] • We use 
this approximation enabling us to compute analytically the liquid density pro- 
files. Nevertheless, we take surface effects and repulsive forces into account by 
adding density functionals at boundary surfaces. In mean-field theory, London 
potentials of liquid-liquid and liquid-solid molecular interactions are 

(fin — — ^ , when r > ai and tp u = oo when r < o\ , 
^is = — J > when r > 5 and ip ls = oo when r < 5 , 

where cu et q s are two positive constants associated with Hamaker constants, 
o~i and a s respectively denote liquid and solid molecular diameters, 5 = \{ cr/+ 
a s ) is the minimal distance between centers of liquid and solid molecules jl]. 

We consider the interaction between a solid surface flat at a molecular scale 
(but curved at several nanometer scale) and a liquid by means of a continuous 
model. The density-functional of energy E of the inhomogeneous liquid in a 
domain D of differentiable boundary S (external forces being neglected) is 
taken in the form 



E 



Ef + Es with E f = J J J pe dv, E s = J J ds. 



The first integral (energy of the volume) is associated with square-gradient 
approximation when we introduce a specific free energy of the fluid at a given 
temperature, e = e(p,/3), as a function of liquid density p and (3 = (gradp) 2 . 
Specific free energy e characterizes together fluid properties of compressibility 
and molecular capillarity of interfaces. In accordance with gas kinetic theory 
[5], scalar A = 2p e^{p, (3) (where denotes the partial derivative with respect 
to (3) is assumed to be constant at a given temperature and 

pe = pa{p) + ^ (grad p) 2 , 

where term (A/2) (grad p) 2 is added to the volume free energy pa(p) of a 
compressible fluid. We denote the pressure term by P{p) = p 2 a'(p). The 
second integral (energy of the surface) is such that the free energy per unit 
surface is [2], 

0(P) = -7iP + \l2p 2 - (!) 
Here p denotes the limit liquid density value at surface S. Constants r y 1 , 7 2 



2 



and A are positive and given by relations [6], 

TTQs TTCu 2-KCu 

71 ~ m 2 m ims Psoh 72 ~ m 2 mf ' 3(7, mf ' 

where m/ et m s respectively denote the masses of liquid and solid molecules; 
p sol is the solid density. 

In this paper, we first develop the boundary conditions for the general case of 
the interaction between a non-homogeneous liquid and a curved solid surface 
with a surface energy due to intermolecular interactions and depending of the 
fluid volume deformation. Then, for a surface energy in form ([T|) we study the 
stress vector distribution on a surface where bumps and hollows are period- 
ically distributed. Finally, we estimate the stress effects for a silicon surface, 
with a curvature of several nanometer range, in contact with water. 



2 Boundary conditions 



The equation of equilibrium and boundary conditions are obtained by using 
the virtual power principle [7118] . For example, virtual displacements £ = 5x are 
defined in a classical way by Serrin [5] page 145, where x = {x 1 }, (i = 1, 2, 3) 
denotes the Euler variables in a Galilean or fixed system of coordinates. 
A liquid (in drop form) occupying a domain D of the physical space lies on a 




Fig. 1. Vector n is the unit normal vector to S exterior to D; vector t is the unit 
tangent vector to T with respect to n; n' = t x n. 

solid surface S (the liquid is also partially bordered by a gas); the edge T (or 
contact line) is the curve common to S and the boundary of D (see Fig. 1). 
All the surfaces and curves are oriented differential manifolds (rj). 



1 Transposed mappings being denoted by T , for any vectors a, b, we write a T b 
for their scalar product (the line vector is multiplied by the column vector) and 
ab T or a <g> b for their tensor product (the column vector is multiplied by the line 
vector) . The image of vector a by a mapping B is denoted by B a. Notation b r B 
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2.1 Variation of the density-functional of energy E 



The density in the fluid has a limit value at the wall S. Then, on S, 

5(f) = (f)'(p) 5p = -p (p'(p) div C, 

where 5p + p div £ = [9] . Let us denote 

G = -p0'(p), H = <f>(p)-p<f>'(p). 

The function H is the Legendre transform of <p with respect to p. For any 
virtual displacement £ null on T, Rel. f[TT|) of Appendix yields, 

where, now S is the imprint of D on the solid surface. Consequently, from the 
calculations in Appendix, we obtain: 

For any virtual displacement null on the complementary boundary of D with 
respect to S and null on the edge T, the variation of E is, 



5E = - J J J (dwa)Cdv 



+ / L.( G - A )^JT + ( o nT + g rad J ( A -H) + n T a)C ds, 



where 



— I — (2) 
ax y ax 

is the symmetric stress tensor of the inhomogeneous liquid, with p = p 2 e, p — 
p div(A gradp) ; A = Xp (dp/dn) with dp/dn = (dp/chc) n ; £ n = n T £ ; 
2/R m is the mean curvature of S and grad tff denotes the tangential part of the 
gradient relatively to S. 



2.2 The virtual work of forces exerted on D 



The virtual work of elastic stresses on S is 

STe = I Is dS ' 

means the covector c T defined by the rule c T = (B T h) T . The divergence of a 
linear transformation B is the covector divi? such that, for any constant vector a, 
(div B) a = div {B a). If / is a real function of x, df /dx is the linear form associated 
with the gradient of / and df /dx 1 = {df /chc), ; consequently, (df/dx.) T = grad/. 
The identity tensor is denoted by I. 
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where k = —a e n is the loading vector associated with stress tensor a e on the 
wall in classical theory of continuum mechanics (rh . Then, the virtual work 
of forces St exerted on D is St = — SE + St c and, 

St = J J J (div<r)Cdv 

ls (G ~A)^+ (^rf 1 n T + gradJ(A - H) + n T a - A C da. 



2.3 Results 



The fundamental lemma of variation calculus, applied to the relation St = 0, 
for all previous virtual displacements, yields: 



The well-known equation of equilibrium for capillary fluids jlOj . 

diva = 0, (3) 
The boundary conditions on S, 

G-A = 0, 
1 n = n + gr&d tg (A - H) + a n. 

Equation (Tj0)i yields a condition relative to the surface energy (CQ) which de- 
pends on the fluid density at the surface and on the quality of the solid wall: 

Xf n+m= , or A^ = 7l - 72 p. (5) 

Equation ([5]) expresses an embedding effect for the liquid density Such a con- 
dition appears for simpler geometry in [1JTT] . 

Condition ^ 2 appears in the literature [7ffT2] but without the terms cor- 
responding to the molecular model ([T]) of surface free energy. Such type of 
condition also appears in interfacial problems with other solid surface energy 
but with a null curvature as in [TyTT] . In Cauchy theory, we are back to the 
classical equation k = a n. 

The definition ([2]) of a implies an = —pn — A (dp/dn) grad p . Then, for an 
elastic wall, by taking into account of Rel. ([5]), the vector k is normal to S, 

rp 1(f) 

k = K n n with K n = n an— — — . (6) 



2 It is important to note that the external unit normal to S with respect to the 
solid is — n. 
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We obtain the stress vector values of the solid at the elastic wall (which is 
opposite to the action r of the liquid on the elastic wall). Relation looks 
like the Laplace formula for fluid interfaces. Nonetheless, we will see in the 
next section some differences between the results for fluid interfaces and for 
liquid-solid interfaces. 



3 An example of elastic effect on a solid surface 

3.1 General considerations 

As bibliography about elastic effects on surface physics, one may refer to the 
review article [T3~| . 

The aim of this section is to present an example of system such that the meso- 
scopic effects of a liquid locally generate important molecular stress vectors 
on a solid surface. We consider a periodic domain such that the substrate 
solid surface has an alternated structure. The solid surface can be considered 
as a flat domain at the Angstrom scale because roughness and undulations 
are only of several nanometer length (such a model is presented on Fig. 2). 
At level with respect to the third axis, the lateral boundary of domain 




Fig. 2. We consider the model consisting of a surface S with bumps and hollows 
periodically distributed on a period L of several nanometers in two directions such 
that, with respect to the third axis, the bump and hollow levels are opposite. Ex- 
trema of the surface mean curvature are located at point A and B; curve C is the 
limit curve of the periodic rectangular parallepiped. Surface S' delimits the liquid 
bulk (at a distance h of a great number of nanometers from surface S). Surface £ 
is the lateral boundary of D. Vector k is normal to S' and z is directed along k. 
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D follows the curve C of the bludging surface. Due to the axial symmetries 
around the lines Ak and Bk, in local coordinates with these lines as third axis, 
gradp = (dp/dz)k and on these lines the stress tensor a of the inhomogeneous 
liquid gets expressions in the form 





d\ 










a = 





a 2 





, with < 










a 3 





ax 
a 3 



«2 



-v 




XpAp 



where A is the Laplace operator. Consequently, on these lines, Eq. 
a constant value for the eigenvalue CI3, 



yields 



p + A 




Pi. 



(7) 



where Pi denotes the uniform pressure in the liquid bulk of density p z bounding 
the liquid layer at level h. 

- Due to symmetries of domain D, we deduce the average stress actions of 
the liquid on S and S' are opposite and numerically equal to the pressure P\. 

- From Rels. (EH) we obtain, at points A and B, a stress vector r = — k, 
action of the liquid on the elastic wall in the same form than the Laplace 
formula form for fluid interfaces, 



Pt + 



Pn 



11. 



- We must emphasize that Rel. (JE1) is only valid at points A and B. In fact, 
Rel. © yields 

2 " • (9) 



n an 




but for points which are not the summits of bumps or the bottoms of hollows, 
— n T crn = p + X (dp/dn) 2 7^ Pe where \{dp/dn) 2 7^ A (gradp) 2 . Consequently 
at a mesoscopic scale, due to the anisotropy of the liquid on curved solid sur- 
faces, Rel. (EI) replaces Laplace's formula of fluid interfaces. 
- The stress vector is directed as k at points A and B. Due to the axial 
symmetries around the surface extrema at points A and B and opposite mean 
curvatures, when we neglect Pi with respect to 2<p/R m , the stress vector as- 
sociated with the hollow corresponding to point A is a vector T parallel to k 
and the stress vector associated with the bump corresponding to point B is a 
vector — T; the two vectors generate a torque on the surface. This result is in 
accordance with results in p3] where the interaction between liquid and solid 
is represented as localized dipoles and monopoles depending on bumps and 
hollows of the surface S. 
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Physical constants 


di 


ui 


mi 


Pi 


Water 


1.4 x 10~ 58 


2.8 x 10 -8 


2.99 x 1(T 23 


0.998 


Physical constants 






m s 


Psol 


Silicon 


1.4 x 10" 58 


2.7 x 10~ 8 


4.65 x 10" 23 


2.33 


Deduced constants 


5 


A 


7i 


72 


Results (water-silicon) 


2.75 x 10 -8 


1.17 x 10~ 5 


81.2 


54.2 



Table 1 

The physical values associated with water and silicon are obtained in references 
[4|15] and expressed in c.g.s. units (centimeter, gramme, second). No information 
is available for water-silicon interactions; we assume that cu = q s . 

3.2 Application to explicit materials 



At 9 = 20° Celsius, we consider water damping a wall in silicon. The experi- 
mental estimates of coefficients defined in Section 1 are presented in Table 1. 
Far from the liquid critical point, the liquid density at the wall is closely the 
same than the liquid density in the bulk [T6] . 

If we consider a mean radius of curvature of surface S, R m = — 10 -6 cm at point 
A and R m = 10~ 6 cm at point B, when we neglect Pi, we immediately obtain 
an arithmetic value of r n = n T r = 10 8 cgs (or 100 atmospheres) corresponding 
to stress effects of large magnitude between areas around points A and B. 

The elastic effects of a liquid on a solid surface result from the topology of the 
contact interface. It is amazing to observe that a solid surface considered as 
an interface between solid and liquid does not require new concept but only a 
supplementary surface energy and likewise surface morphology. 
An important assumption in the previous calculations is that three scales 
infer in the surface physics: a length scale of one nanometer associated with 
molecular effects and the expression of surface energy, a length scale of ten 
nanometers associated with the size of undulations and surface roughness and 
a length scale of one hundred nanometers associated with the distance of the 
liquid bulk to the surface S. 
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4 Appendix 



Let S be a differentiable oriented manifold in the 3-dimensional space and n 
its oriented unit normal locally extended in the vicinity of S by the expression 
n(x) = grad d(x), where d(x) is the distance of point x to S; covectors grad T a 
and grad^a denote the transposition of grad a and grad tff a, respectively; for 
any vector field w, we get [7j: 
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/xt t 9n dw 

rotfn x w) = ndiv w — wdivn + — — w — — n. 

ax ax 
dn 2 

From n T — — = and divn = we obtain on S: 

ax R m 

T I \ 1 • 1 rp rp 9"W 

n rotfn x w) = div w H n w — n — — n, 

Rm OX 

and we deduce, 

Lemma 1 For any differentiable scalar field a, 



(10) 



where gracL a 



a div £ = a 

da 



jjT (n ~ (gradja) C + n T rot (an x £). 



dn <9x 



<9x 



/ — nn J 



belongs to the cotangent plane to S and 



Application to the calculation of 8E 



All the densities are expressed in the physical space. The domain D is a ma- 
terial volume [9], then bEf = J J J p Se dv. 

_ . de de . _ . . dp db p dp d£ , rTTTIN 

From be = — bp + — bp and b — = — — — (see pjJJ), we get: 

op dp ax ax ax ax 



pe )/3 8(3 = 2pe,p5 



dp I dp s 
<9x \ <9x 



A 



' dbp dp dC,\ [dp' 
<9x 9x9x/ \ <9x 



T 



= div(A grad p bp) — div(A grad p) Sp — tr I A grad p grad p — — 



Due to o~p = — pdiv£ (see [9]), with a = Xp(dp/dn) = A and by usinj; 
definition O), we obtain: 



/ dp 

D \ OX 



div (A grad p grad T p) 1 C,dv 



div ( A p grad p div £ + A grad p grad T p C + P C ) dv 



ds. 



= J J J D -(diva)Cdv + J J {-A divC + n^C) 

From the Stokes formula, we get: 

n T rot (An x Q ds = [ A t T (n x O di = [ A n' T C d£, 
Jr Jr 
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which is null in the case of the virtual displacements of Section 2.1. Finally, 
by using Rel. (ITTi) . 

"W / I -< d - « **/ /, { f + (f nT + grad s 4 + nIj ) <} di 

• Application to the calculation of 5E S : 



Due to Eg = J J det (n, d]X, d 2 x) where d]X and <i 2 x are two coordinate 
lines of S, we get: 



E s = / / det F det (F _1 n, diX, d 2 X), 

where Sq is the image of S in a reference space with Lagrangian coordinates 

X and F is the deformation gradient tensor — — of components < — — r > , (see 

<9X 1 9X3 J 

|S]). Then, 

<5F S = y jf 50 detFdet(F~ 1 n,d 1 X,d 2 X)+y > jf 5 (det F det(F _1 n, diX, d; 
with [ i (f>6 (det F det (F _1 n, diX, d 2 X 

9n \ / <9£ 



X~ 



So 



J J (ft div C det(n, dix, <i 2 x) + 4> det ^— d]X, d 2 xj — det ^— n, d]X, d 2 x 
| jf (div(0 C) - (grad T 0) C - n T || ds. 



Relation (PSJ yields: 



2 _ T/> _t^0 C r 



div (0 C) H n (-n — — n = n rot (0 n x £). 

Rm OX 

Consequently, 

• 5 (det FdetfF-^diX^aX' 

. v ' 

^ n T + grad T (nn T — I) J £ ds + / / n J rot (0n x £) ds, 



's'o \ F m y J Js 

and finally due to ^ n /T £ d£ = for the virtual displacements of Section 2.1, 



2: 



<5F S = / / (50-(^n T + gradj0 C ) ds. 
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